Emission properties of a driven artificial atom: increased coherent scattering and 

off-resonant sideband narrowing 
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We study the first order field correlation function and emission spectrum of a driven artificial 
atom. Focusing on the experimentally relevant example of a semiconductor quantum dot, we find 
that the solid-state (phonon) environment can have a profound effect on the emission characteristics 
of such a system. For resonant driving, we find an increase in the coherently emitted radiation field 
with increasing driving strength due to the quantum nature of the phonon bath. This feature is in 
stark contrast to the standard quantum optics treatment, which predicts a decrease in the fraction 
of coherent emission with increasing driving strength, and should be observable in experimentally 
achievable regimes. For off-resonant driving, we show how narrowing in the spectral sideband width 
can occur in certain regimes, consistent with a known experimental trend. 



As described by Mollow, the temporal correlations and 
fluorescence spectrum of light scattered from a driven 
two-level system depend crucially on the relative size of 
the driving strength and the two-level system radiative 
decay rate [l|. At weak driving strengths, the light is 
predominately coherently (or elastically) scattered, re- 
sulting in the absence of oscillations in the first-order 
field correlation function, and a single peak in the emis- 
sion spectrum at the driving laser frequency. At larger 
driving strengths, however, the coherently scattered light 
is strongly suppressed, and the emission becomes dom- 
inated by incoherently scattered light. This results in 
oscillations in the field correlation function and a triple- 
peak structure in the spectrum, known as the Mollow 
triplet. 

While these fundamental predictions have long been 
confirmed in the traditional quantum optical setting of 
driven atoms Q, they have only much more recently 
been observed in solid-state artificial atoms such as 
driven semiconductor quantum dots (QDs) yHSII , single 
molecules [iflj, and superconducting circuits In fact, 
in the particular case of QDs, many of the archetypal fea- 
tures of atomic quantum optics have now been seen, such 
as resonance fluorescence [3j-|9|, coherent population os- 
cillations 0, [H-14 1, photon anti-bunching 15, 3], and 
two-photon interference Aside from being of fun- 

damental interest, these observations constitute signifi- 
cant progress towards the ultimate goal that QDs may 
be used as efficient single photon sources US-SI, and for 
other quantum technologies [241 ]. 

Though a driven QD thus bears a great deal of resem- 
blance to the more idealised case of a driven atom in 
free space, it is nevertheless unavoidably coupled to its 
surrounding solid-state environment. For ground-state 
excitonic transitions in typical InAs QDs, coupling to 
acoustic phonons has been demonstrated to dominate the 
QD-environment interaction 13,3]. This causes the be- 
haviour of the QD, and hence also of its emission prop- 
erties, to deviate from the atomic case. In particular, a 



now well-known signature of phonon coupling is the ap- 
pearance of an excitation-induced dephasing contribution 
which depends on the square of the Rabi frequency, with 
a rate 7 ~ fl 2 [1, [H , 14 , 25| . This driving dependence is 
theoretically understood as resulting from phonons that 
induce transitions between the dressed states of the QD 
at the Rabi energy 2|§- 28| , making it the relevant energy 
scale in the three-dimensional phonon environment. 

Thus, to accurately describe the emission properties 
of a driven QD, one must faithfully include its coupling 
to phonons. This can be achieved in various ways, such 
as througj^ master equations of weak-coupling [25|, |26| |. 

*32j | . as well as by 
For our purposes, 



polaron 28-31], and variational tyj 
several numerical methods 27. 33l 3 
master equation techniques are particularly attractive 
since, with use of the quantum regression theorem 35|, 
they can readily be extended to investigate field correla- 
tion properties. For example, in Refs. 3^, 3l| a polaron 
master equation was used to explore the incoherent spec- 
trum of a QD in a micro-pillar cavity, while QD emission 
under both pulsed [36 , 37 1 and incoherent [38| excitation 
has also recently been studied. 

Here, we extend the variational method of Ref. [32[ to 
explicitly investigate phonon effects in both the coherent 
and incoherent emission properties of a QD, driven on 
and off resonance. On resonance - and in stark contrast to 
the atomic case - we find that at moderate driving the co- 
herent contribution to the QD resonance fluorescence can 
actually increase with driving strength. This results from 
an imbalance of phonon absorption and emission pro- 
cesses at low temperatures, an effect that arises naturally 
in our miscroscopic model of the phonon bath, but can- 
not be captured by a simpler treatment in terms of a phe- 
nomenological pure dephasing process. We also find that, 
in an appropriate parameter regime, our model predicts 
a narrowing of the Mollow sidebands as the QD-laser de- 
tuning is increased, which is consistent with recent (and 
as yet unexplained) experimental observations [8J. 

Methods summary - We model the QD as a two-level 
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system with ground state |0) and excited (single exciton) 
state \X), split by an energy haj . The dot is coherently 
driven by a laser of frequency w/, with Rabi frequency Q, 
and is coupled to two separate harmonic oscillator baths 
to account for both phonon interactions and spontaneous 
emission into the radiation field. After moving into a 
frame rotating at frequency cj;, and making a rotating 
wave approximation on the driving term, our Hamilto- 
nian takes the form (for h = 1) 

H = v + ~~(t x + wk&£&k + wl** 

k q 

+ \X)(X\ ]T g k (bt + 6 k ) + ^(/ lq a q c^ t ( 7 + + H.c), 

k q 

were v = cj — is the QD-laser detuning, cr + = |AT)(0| 
(cr_ = <r+), g x = cr + + <7_, and H.c. stands for the Her- 
mitian conjugate. The phonon bath is represented by 
the creation (annihilation) operators (&k) for modes 
with frequency wt, and couples to the QD through the 
constants <?k- The photon bath is similarly defined, with 
operators a q (a q ), frequencies 7y q , and coupling constants 
hq. 

Obtaining an equation of motion for the QD degrees of 
freedom can be achieved by several means. As mentioned 
previously, in this work we opt to use the variational mas- 
ter equation method of Ref. [32J , since it has been shown 
to be unrestricted to either weak phonon coupling, or to 
the small driving limit of polaron theory (2|| ; it therefore 
constitutes a robust yet efficient way in which to include 
the effects of phonons within the QD dynamics. In brief, 
the variational master equation is derived by first apply- 
ing a unitary (QD state-dependent) displacement trans- 
formation to the full QD-phonon Hamiltonian, which is 
subsequently optimised on the basis of free-energy min- 
imisation arguments |39| . The second-order, time- local 
projection operator technique is then used to approxi- 
mate the QD dynamics in the transformed representa- 
tion. For our present purposes, the great advantage of 
using a combination of the quantum regression theorem 
and a master equation obtained in this way, is that the 
Born approximation on which the regression theorem re- 
lies [35j is performed in a transformed frame, and is not 
therefore restricted to weak phonon coupling. Full de- 
tails of the derivation - extended in this work to include 
the photon bath - can be found in the Appendix. 

We find that the Schrodinger picture master equation 
takes the form 



is the phonon-renormalised 
with P = l/fc^T the invers 



Rabi frequency, 
temperature, and 



Pv 



[ea z + Q, r a x ,p v ] +/C p h(/3v) + £ sp (py), (1) 



where py is the reduced density operator for 
the QD two-level-system in the transformed 
frame, e is the detuning of the laser from the 
variationally shifted QD transition frequency, 
n r = ttexp[-i r o °° J ph (w)uj- 2 F(uj) 2 coth{f3uj/2)dLu} 



fc sp {pv) — r , i(<r_/0y<7 + — (l/2){er + cr_, p v }) accounts for 
spontaneous emission processes at a rate T\. The func- 
tion Jph(w) = J2k dk^i^^^k) is the QD-phonon spectral 
density, which for coupling to acoustic phonons is ade- 
quately described by J p h(w) = auj 3 exp[— (uj/uj r ) 2 ], with 
coupling constant a and cut-off frequency oj c [13l Il4 |. 
The factor F(uj) arises from the variational transfor- 
mation used in deriving Eq. (fT]), becoming unity in 
the full polaron limit. The term /C p h(/?v) contains all 
phonon effects other than those included in the renor- 
malised Rabi frequency fl r , and in general has a fairly 
complicated form, representing the various processes 
induced by phonon interactions (e.g. pure-dephasing 
and phonon emission/absorption). We shall see below, 
however, that in certain regimes its dominant effects can 
be approximated by simple intuitive forms. 

Though Eq. (fT]) describes the QD excitonic degrees 
of freedom, by relating the field operators to those of 
the QD it can be used to calculate the steady-state 
first order field correlation function of the QD emis- 
sion, <7^(t) = lim t _ i . 00 (a + (t)a-(t + r)) (35|. Impor- 
tantly, g^(r) does not necessarily approach zero as 
t — > oo, and in this limit gives the coherent contri- 
bution: g^jj = linv^oo g^ (r). The coherent contribu- 
tion can be shown to be equal to the square modu- 
lus of the off-diagonal elements of the QD density op- 
erator in the steady-state, g^ h = \pox\ 2 , and is thus a 
direct consequence of non-vanishing coherences in the 
QD. The difference between the full correlation function 
and the coherent part is then the incoherent component, 



5inc( r ) = 9^K T ) ~ 5coh- Armed with the first order field 
correlation function, we may then calculate the QD emis- 
sion spectrum, S(w) cx (l/?r)Re[ / °° e l ^~^ T g^ {t)At]. 

Enhanced coherent scattering - We begin our analysis 
by investigating the emission properties of the QD when 
driven on resonance with the polaron shifted transition 
frequency [8 = v — J2k dt^k 1 = 0)- We are interested in 
examining the detailed effects induced by the coupling to 
phonons as the driving strength is varied. In particular, 
we would like to explore deviations from the more phe- 
nomenological treatment of environmental interactions 
(beyond radiative decay) as giving rise simply to sources 
of pure dephasing j35|. As mentioned, the particular form 
of /Cph(pv) already tells us that in general phonon cou- 
pling should have an effect on the QD emission proper- 
ties beyond the pure dephasing model. However, it is 
instructive to first consider to what extent these effects 
can be captured by simpler terms in the master equa- 
tion. For example, at weak driving strengths (/30 <§; 1) 
we find that we may indeed approximate the phonon 
influence by a pure dephasing form (see Appendix for 
details), /C ph (pv) ~ (1/2)7pd(o- z Pvct z - pv) in Eq. ([T]), 
with a rate given by the expression from polaron the- 
ory, 7p D = (OV 2 ) 2 fZc C0S ( n rs) (e^-e-^ds, where 
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2 (cos(ws) coth(/3w/2) — is'm(ujs))duj, 
2S|, |30(. Thus, at sufficiently weak 



and F(lj) — > 1 in fi 
driving strengths, we find that phonon interactions can 
be represented simply as fluctuations in the excited state 
energy of the QD, and as such give rise only to an ad- 
ditional pure-dephasing process, consistent with experi- 
mental results in this regime @, @, i, [U, [13] ■ With the 
further (minor) approximation that e = in Eq. ([T]), we 
then arrive at an analytic expression for the first-order 
field correlation function of a QD interacting with acous- 
tic phonons, in the weak-driving limit: 



r (ns) 



2fi2 + 21^2 

+ c - i Ti +r 2 )r ( C0S (£ T ) _ M sin(Cr)) 



(2) 



where T 2 = ±r 1 + 7 p D) C = V^r - WW~i ~ r 2) 2 , ^ = 

(n? - r x (Ti - r 2 ))/(2fi= + 2r 1 r 2 ), and m - (^(r 2 - 

3ri)+r3r2(rr 1 -r 2 - 1 ) 2 )/(4C(^ + rir 2 )). The coherent 
contribution in this limit is given by 



o (1) 

i/coh 



2T 1 T 2 + 2fl 2 r 



(3) 



which would clearly approach zero if Q r were allowed to 
become large. 

Eqs. @ and ^ are essentially the standard atomic 
gWfY) expressions, when extended to include pure- 
dephasing 0, E|. The only difference here is that we 
explicitly include a driving dependent pure-dephasing 
rate through 7pd, and that the driving is itself renor- 
malised by the phonon environment through il r . Both 
of these features are important for these expressions to 
adequately approximate the full dynamics. For weak 
exciton-phonon coupling, we can expand the exponen- 
tial factors in the pure-dephasing rate to further ap- 



proximate as 7pd 
fiCl < 1 and n r /a>, 



7pd ~ [13|, 111 UM 



« f J(n r ) coth(j8fi r /2), which for 
<C 1 takes on the well known form 
28] . In the frequency domain, this 
form of pure-dephasing will be apparent in the broad- 
ening of the spectral sidebands with increasing driving 
strength [1, El ED] . 

Thus, at weak driving, phonon coupling does not 
induce behaviour markedly different from the atomic 
case. However, to exemplify the breakdown of the pure- 
dephasing model at stronger driving strengths, in Fig. Q] 
we plot g^'{r) calculated using the full variational the- 
ory (solid blue curves) and calculated using Eqs. ([2j and 
([3]) (black dashed curves). As expected, for low driving 
strengths, f2 < O.lps -1 (< 15 GHz in linear frequency) 
the pure-dephasing model gives an excellent approxima- 
tion to the full theory. As the driving strength is in- 
creased, however, significant discrepancies soon begin to 
become apparent. In particular, from the different long- 
time values that are approached when Q > 0.33 ps" 1 , 
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FIG. 1: First order field correlation function for various driv- 
ing strengths as indicated, calculated using the full variational 
theory (blue solid curves), and calculated using the pure- 
dephasing approximation of Eqs. @ and © (black dashed 
curves). The right-most parts of the upper three plots show 
enlargements of the long-time behaviour. The lower plots 
show only the coherent contribution, now as a function of 
driving strength, calculated using the same methods as in the 
other plots, and also from the simple expression of Eq. Q (or- 
ange dotted curve) . The lower left plot shows the only region 
where the pure-dephasing model gives a non-negligible coher- 
ent contribution, close to the origin. Parameters: Ti = 700 ps, 
a = 0.027 ps 2 , uj c = 2.2 ps -1 , and T = 4 K. 



we can conclude that the coherent contribution begins 
to become important, and that this is not captured 
by the pure-dephasing approximation. Indeed, when 
fl = 4 ps" 1 we see that in the full phonon theory a 
long-time limit of gj^ ~ 0.25 is approached, in clear 
distinction to the pure-dephasing model. 

That Eqs. ^ and §5§ cannot capture these effects sig- 
nifies that above a driving strength of fi ~ 0.1 ps -1 (for 
these realistic parameters), the field correlation proper- 
ties of the QD emission fundamentally depart from the 
atomic case. This departure can be attributed to phonon- 
induced coherences in the steady-state of the QD. The 
pure dephasing model is unable to capture this effect, 
as it does not lead to the correct equilibration with re- 
spect to the phonon bath; it essentially assumes an in- 
finite temperature limit, and thus the quantum nature 
of the environment is lost. However, for resonant driv- 
ing, we may still approximate this process in a sim- 
ple way, by modifying the phonon term in Eq. (JT|) to 
£ P h(pv) ~ (l/2)7PD(e^pvCT z -pv) + (?7 4 ) K K> {<Tz,Pv}], 
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where k = {n r /2) 2 J^ oc sm(n r s)(e^ - e'^ds, such 
that k/7pd = tanh(/3f2 r /2) 28]. With this adjustment 
to gW(r) depart from the atomic 



both contributions 
case 



41]. Specifically, the coherent contribution now be- 



comes 



r 0) - 

^coh — 



r x n r 



2r!r 2 + 2^ 



Bn 



ri + 27p D 



(4) 



where B — Q r /Cl is the renormalisation factor, arising 
here since we must transform out of the variational frame 
to obtain the QD coherences. 

We see that the coherent contribution now has two 
terms. The first, which was present in the pure-dephasing 
model, has a maximum when Q r = y/TiTsf but quickly 
tends to zero as f2 r is increased. The second term is 
present nowhere in the standard quantum optical treat- 
ment and becomes important for larger values of Q r . In 
fact, once the driving is strong enough such that T± <C 
7pd, we can approximate G^ h s» {B tanh(/3f2 r /2)/2) 2 . 
Thus, with reference to Fig. [TJ we see that at a tem- 
perature of 4 K, we have ksT s» 0.5 ps , and when 
driving at fi — 0.1 ps -1 we are still effectively in a high 
temperature regime (/3f2 < 1). As such, the coherent con- 
tribution arising from the QD-phonon coupling is unim- 
portant, and the pure-dephasing approximation is ade- 
quate. As the driving is increased, however, the quan- 
tum nature of the phonon environment begins to play 
an increasingly important role, and the coherent contri- 
bution is enhanced correspondingly. When the driving 
becomes very large, the coherent contribution saturates 
at around B 2 /4 |42j. These observations are borne out in 
the lower part of Fig. [TJ where the coherent contribution 
is plotted as a function of Cl using the full variational the- 
ory together with the approximate expression of Eq. (j4j) , 
which captures the general trend, though gives a slight 
underestimate at large driving. The lower left plot shows 
the region close to the origin, the only regime in which 
the pure-dephasing model gives a non-negligible coherent 
contribution. 

Spectrum - Having investigated the behaviour of the 
coherent contribution to the QD emission, we now turn 
our attention to the properties of the QD emission spec- 
trum. We shall now concentrate on the regime where the 
incoherent contribution dominates (i.e. relatively weak 
driving), and Eq. ([2]) is thus approximately valid on res- 
onance. From a Fourier transform of Eq. ([2]), we find that 
the resonant Mollow sideband widths are determined by 
Ti + r 2 = (3/2)Ti + 7pd, while their positions are given 
approximately by ±Q r . As such, we expect them to show 
a systematic broadening and splitting with increasing 
driving strength, as has already been confirmed [j, l30l]. 
Off resonance, we might then also expect sideband broad- 
ening and splitting with increasing detuning (for fixed f2) 
if we were to replace fl r with the generalised Rabi fre- 
quency, Q' r — y / T^r+~5 2 @- This would lead to similar 
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FIG. 2: From top to bottom, incoherent emission spectrum, 
extracted sideband splitting, and extracted sideband width 
for varying driving strengths on resonance (left), and varying 
detunings (right). The solid black curves in the emission spec- 
tra are for 5 — 0, and a driving strength of f2 = 0.025 ps - 
(which sets our x-axis units in the rest of the plots). The 
dashed red curves are for — 0.094 ps -1 on resonance, and 
S = O = 0.025 ps -1 off resonance (which has been enhanced 
by a factor of 5). The insets in each case show the red side- 
bands shifted and rescaled to lie on top of each other. The 
solid blue curves in the middle row show the functions 2Q r 
(left) and 2^/0.'^. + S' 2 (right). The symbols in the bottom 
row correspond to the red (x) and blue (+) sidebands. Pa- 
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T =10 K. 
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trends for increasing 8 as for ft. However, the experi- 
ments of Ref. Q showed a systematic narrowing of the 
Mollow sidebands with increasing detuning, despite the 
fact that the splitting is well captured by a generalised 
Rabi frequency of a form similar to Vl' r , leaving open the 
question as to why this might be the case. 

In fact, off-resonance, the expressions for the spectrum 
become significantly more complicated than in the reso- 
nant case, and the above simple reasoning does not hold. 
To illustrate this, in Fig. [2j from top to bottom, we plot 
the incoherent emission spectrum, extracted sideband 
splitting, and extracted full-width-half-maximum (= T) 
of the Mollow sidebands, calculated from the full phonon 
theory. In the latter two cases, the spectrum was fitted 
to a sum of three model Lorentzian functions of the form 
L(w) = 0.5r/[(w - uj v f + (0.5r) 2 ]. The left column cor- 
responds to varying the driving frequency on resonance, 
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while the right column corresponds to varying the detun- 
ing with a fixed driving strength. For the spectra in the 
top row, the black solid curves both correspond to driv- 
ing at a strength = 0.025 ps _1 on resonance, while the 
red dashed curve on the left corresponds to increasing the 
driving strength up to £1 = 0.094 ps" 1 , and that on the 
right has values of S = O = 0.025 ps _1 . Before any de- 
tailed analysis is performed, from the shape of these two 
emission spectra, it is evident that moving off-resonance 
is not equivalent to increasing the driving strength; the 
asymmetry in the off-resonant case cannot be captured 
by any combination of parameters in the resonant case. 

A more detailed analysis of the spectra reveals two 
important trends. Firstly, as can be seen by the side- 
band splittings in the middle row, increasing the driving 
strength on resonance does, as expected, cause the side- 
bands to move apart linearly with f2. Also, we see that 
moving off-resonance appears to alter the sideband split- 
ting in exact accordance with the simple procedure of re- 
placing O r — > \J + S 2 . The extracted sideband widths 
in the lower row, however, reveal something quite differ- 
ent. On the left we plot the extracted sideband widths 
as a function of driving strength on resonance. In accor- 
dance with our approximate pure-dephasing model, for 
which we have 7pd ~ , we see a systematic broadening 
of the sidebands with increasing driving strength. On the 
lower right, we plot the extracted widths of the red (red 
crosses) and blue (blue crosses) sidebands as a function of 
detuning. In contrast, we now see a systematic narrowing 
of the sidebands [43} , consistent with recent experimental 
results Q. To further confirm this point, the insets of the 
plots in the top row show the red sidebands in each case 
plotted on top of each other. It is important to note that, 
while we have attempted to model neither the effects of 
the cavity in Ref. which can be significant [3^, 31 1, 



nor other sources of noise beyond phonons, these results 
still demonstrate that an increase in sideband splitting 
off-resonance does not necessarily lead to an associated 
phonon-induced increase in sideband width. 

We can gain some approximate analytical insight into 
this behaviour by again considering the pure-dephasing 
limit. Allowing for off-resonant driving, we expand the 
sideband widths to second order in the detuning, from 
which we find that they are now determined by (3/2)ri + 
7pd - {5/\/2Q. r ) 2 (Ti - 27pd). Hence, for Ti > 27 PD , as 
is the case in Fig. [2j we expect narrowing with increas- 
ing detuning from resonance. Finally, it is also worth 

noting that the pure dephasing rate is itself modified off- 
si 2 

resonance, and becomes 7pd ~ § jpr J(fl' r ) c °th (/3Q' r /2) 
in the weak-coupling limit. For our cubic spectral den- 
sity, we then find 7pd ~ aCl^Q' r exp [— W?/cj%\ for low 
temperatures (j3fl' r ^> 1) which increases with detun- 
ing for Vl' r -C uj c . On the other hand, at high tempera- 
tures {fiil' r < 1) we have 7PD - a(fl 2 r /(3) exp [-fl' 2 /oj 2 }, 
which decreases with increasing detuning, though slowly 



for n' r <C lo c . This is the relevant regime for Fig. [2] 

In summary, we have studied the coherent and incoher- 
ent contributions to the first order field correlation func- 
tion and emission spectrum of a driven artificial atom, 
focusing on the role played by the solid-state (phonon) 
environment. For resonant driving, we find an increase 
in the coherently emitted radiation field with increas- 
ing driving strength due entirely to the influence of the 
phonon bath. This feature is thus absent in the standard 
quantum optics treatment, which predicts negligible co- 
herently emitted radiation in the same limit. For weak 
off-resonant driving, we have also shown that despite the 
fact that the generalised Rabi frequency increases with 
detuning, phonon interactions do not necessarily lead to 
a related increase in spectral sideband width. In fact, 
sideband narrowing can occur in certain regimes, consis- 
tent with an observed experimental trend [8] . Though we 
have concentrated here on the case of a driven QD, the 
model we have employed is sufficiently general to apply 
to a wide range of two-level systems in which environ- 
mental interactions (beyond spontaneous emission) are 
significant. We would therefore expect the qualitative 
features of artificial atom emission that we have outlined 
above to be observable in various experimental settings. 
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APPENDIX 

In this Appendix we give the derivation of the master equation used in the main text, Eq. (1). We first show how 
the quantum dot-phonon and quantum dot-photon effects can be treated independently within our formalism. We 
then give the resulting expressions of the variational method used to treat the quantum dot-phonon coupling, and 
shown how they reduce to a pure-dephasing form in the appropriate limit. 

The system 

We model the quantum dot (QD) as a two-level system with ground-state |0) and excited (single-exciton) state 
\X), separated by an energy Hujq. The QD is continuously driven by a laser of frequency u>i, with Rabi frequency f2, 
and is coupled both to a phonon bath, represented by an infinite collection of harmonic oscillators of frequencies 
and creation (annihilation) operators 6 k (6k), arid to the radiation field, again represented by an infinite collection 
of harmonic oscillators, this time of frequencies ry q and creation (annihilation) operators a q (a q ). The two baths are 
assumed to be entirely independent, such that [a q , bjj = [6k, a q ] = 0, etc. 

In a frame rotating at the laser frequency, and within the dipole and rotating wave approximation for the QD- 
radiation field coupling, the complete Hamiltonian has the form (setting h = 1) 

H = u\X){X\ + ^a x + £> k & k &k + \X)(X\ ^9^1 + &k) + J2 ^ a q a q + ^(Mq^+e^* + h^a.e-^), (5) 

k k q q 

where a rotating-wave approximation has also been made on the driving term. Here, v = luq — u)\ is the detuning of 
the laser from the excitonic transition energy, the QD-phonon couplings are denoted by c/k (assumed real) and the 
QD-radiation field couplings by /i q . The Pauli matrices are defined in the usual way, such that a x = <r_ + er + = 
|0)(A| + \X)(0\ and a z = \X)(X\ - |0)(0|. 

Variational transformation 

We now apply a unitary variational transformation to H in order to attempt to treat the QD-phonon interaction 
beyond the weak coupling approximation [32j]. The transformed Hamiltonian is defined by Hy = e v He~ v , where 

V=\X)(X\J2(^bt-aib k ), (6) 

k 

with Ok = /k/wk, and where /k are variational parameters to be determined later. We can write e ±v = |0)(0| + 
| X)(X\ Y[ k D(±ak) where f| k £'(±ak) = Ilk ex P[ = ' = ( ak ^k — a k^k)] is a product of displacement operators D(±ak). 
After the variational transformation, we find Hy = H$ + H\ a + Hu, + Hyi + Hb 1 + Hb 2 , where 

H S = + -a z + —a x , (7) 

the interaction terms H\ a — \X){X\ ^2 k (gk — fujib^ + 6k) and Hu, = ^ {&xB x + <JyB y ) contain only QD and phonon 
operators, the interaction term 

H 12 = ^(fe q a q B+^e m,t + h^B^a-e-^), (8) 
q 

contains QD, phonon, and photon operators, and the bath Hamiltonians are Hb x — Xk^k^k^k and Hb 2 = 
S q r lci a \ a ci- The detuning is now e = oj' q — wj, defined in terms of the bath-shifted QD transition energy lu' q = ojq + R, 
with R = X)k w k VkC/k — 2t?k)- The bath operators are given by B x — \{B + + B- — 2B) and B y — -^(B- — £>+), 
where B± — J^ k Z?(±ak), which have the same average with respect to a thermal equilibrium state of the phonon 
bath: B = (B±) = cxp[— (1/2) ^ k |ak| 2 coth (/3wk/2)], with inverse temperature (3 = l/(/csT). The bath-renormalised 
Rabi frequency, which appears in the system-Hamiltonian as a result of our choice of bath operators B x and B yi is 

r^ r = b^i. 
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Time-local master equation 



To proceed with the master equation derivation, let us now separate the variationally-transformed Hamiltonian into 
interacting and non-interacting terms such that Hy — Hq+Hi, with Hq = Hs + Hb 1 + Hb 2 and Hi = Hi a + Hi b +Hi2, 
and treat Hi as a perturbation. Notice that while the interaction terms H\ a and -Hit, depend only on QD and phonon 
operators, and not on the radiation field, the term H\2 does not depend solely on QD and radiation field operators, 
but also on the phonon operators B± as a result of the variational transformation. At this stage, we could perform a 
mean- field- type approximation and replace B± by their average B in Eq. ([8"j), which would significantly simplify the 
form of H\i- However, we shall not perform such an approximation here, but instead attempt to account explicitly 
for the interplay between the phonon and radiation field operators in the following master equation derivation. 

We now move into the interaction picture with respect to H$, yielding an interaction Hamiltonian in the 
(variationally-transformed) interaction picture of the form Hx(t) = Hi a (t) + H\ b + Hi 2 (t), where H\ a {t) — 
e lHot H la e- lHa \ H lb (t) = e iHot H lb e- lH °\ and H 12 (t) = e lHot H 12 e- lHot (tildes indicate interaction picture op- 
erators). We find H la {t) = 7fx(*)Ek(fk- /k)(6] c e ltJk * + b k e~ iUkt ) with n x (t) = e iHst \X){X\ e~ iHst , while 
H lb (t) = S (a x {t)B x (t) + a y (t)B y (t)) where a n (t) = e tHst a n e- tHst and B n (t) = e iHB ^B n e- iHB ^ for n = x,y. 
Explicitly, we may write B x (t) = \{B+{t) + B_(t) - 2B) and B y {t) = ^(B-(t) - B + (t)), where B±{t) = 
e iH Bl t B±e -iH Bl t _ ]T k £)(±ake <Wkt ). The third term (QD-photon-phonon) is given by 

H 12 (t) = ^(/ lq a q e-^* J B + (t)<T + (t)e^' t + h^e^B^^^e'^), (9) 
q 

where a±(t) = e lHst a±e~ lHst , and B±(t) are as defined above. We now consider the products 



ff±(0e ±iUlt = exp [i(~a z + ^-a x )t]a ± exp [ - i(~* z + y**)*] e ±i{u °- v)t (10) 

which appear in Eq. ([9]). Provided that ujq 3> u, e, Q r , which is generally the case for driven QDs since ujq ~ 1 cV 
compared to meV or smaller energy scales for the other quantities, then we may approximate a±(t)e ±lu '' t ~ cr±e ±l " ot , 
which simplifies the third interaction term to 

H 12 {t) a ^(/i q a q e-^*i3 + (i)a + e^ t + h^aX^B-^a-e^ *), (11) 
q 

without any significant loss of accuracy. 

We now follow the standard projection-operator procedure (ioj to derive a time-local master equation for the 
reduced QD exciton density operator, pv, in the variational frame interaction picture, treating Hi{t) to second order 
(the first order term in Hi{i) disappears). In doing so we arrive at the usual expression 



dpy(t) 
dt 



= - f d S ti B [Hi(t),[Hi( S ),p v (t)p B }}. (12) 
Jo 



Here, p B is a reference state of the environment (comprising both phonons and the radiation field), which we take 
to be a thermal equilibrium state: p B = e~^ B /tr B (e~' 9HB ) = [e~^ B i /tr Bl (e^^i )] [e~P HB 2 /tr Sa ( e -P HB 2 )] = 
Pb 1 Pb 2 - Inserting Hi(t) = H\ a (t) + H\b{t) + ffia(t) into Eq. (fT2|) we find that the master equation consists of nine 
terms corresponding to the different possible combinations of the three interaction terms. However, since {H\ a )B 1 = 
{HibjBt = {H\2)b 2 = Oj we fi n d that those terms containing H\ a or Hu and Hi 2 disappear. We are then left with 



t rt 
dstr Bl [H la (t) + H lb (t),[H la (s) + H lb (s),pv{t)p Bl }} - / dstr Bl+B2 [H 12 (t) ,[H 12 (s), p v (t)p Bl pB 2 ]]- 



dpy(t) _ 

dt 

(13) 

The first term on the right-hand of Eq. (fT3|) is precisely that expected from the variational treatment of QD exciton- 
phonon interactions in the absence of the radiation field, and is thus unaffected by the extra coupling to the elec- 
tromagnetic field. However, the second term, responsible for photon emission and absorption processes due to the 
radiation field, appears at this stage to be modified by the phonon environment, due to our use of the variational 
transformation. We shall now show that for typical QD parameters this term should actually reduce to the usual 
form expected in the absence of phonon interactions. 
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Spontaneous emission terms 



We now write ^gff = Kph(pv(t)) + fc S p{p~v(t)) where 



dt 

ft 



jC sp (p v (t)) = - [ dstv Bl+B2 [H 12 (t),[Hi2(s),psp{t)p Bl PB 2 }}. (14) 
Jo 

and consider this term in more detail. The QD-phonon coupling term ]C p h(pv(t)) is discussed later. To proceed, 
we write H 12 {t) = A{t)Q{t)B + {t) + {A(t)Q(t)B+{t))\ with A(t) = cr+e 1 ^' 1 and Q(t) = E q /l q a q e_ "' q *: where ^( f ); 
Q(t), and B + {t) all commute. Inserting this into Eq. (fl4l) and taking the trace over the environmental degrees of 
freedom, we find that the radiation field term can be written in the simple and familiar form (still in the variational 
representation interaction picture) 

K Bp (pv(t)) =T(t) (a-p v {t)a+-^{a + <T_,pv{t)}\ . (15) 

where we ignore absorption and stimulated emission processes under the assumption that no thermal photons exist 
at the appropriate energy scale for temperatures of interest. Additionally, we have ignored the Lamb-shift of the 
excitonic energy splitting induced by the radiation field. Importantly, the rate of spontaneous emission processes is 
given here by 



t 



r(t) = 2Re/ dre 4 " 0T C(r) V l/iqpe""'^, (16) 
where 

r r°° J , 1 

(17) 



C(t) = exp 



dw Jph ^ F H 2 ((l - coswt) coth/Jw/2 + ismuir) 



o 



with /k = gk-F(wk), is the usual polaron-type bath correlation function for the phonon operators, i.e. C(t) = 
(B±(t)B^(0)) Bl , written in the continuum limit, with J p h(w) = ^ k |<?k| 2 <5(o; — Wk) being the spectral density relevant 
to the phonon bath. Taking the continuum limit for the radiation field in Eq. ()16[) . we can then write the spontaneous 
emission rate as 

r(i) = 2Re f dTe luJ0T C(T)X(T), (18) 
Jo 

where we define the radiation field correlation function X(r) = J °° drye _I,?T J P t(?7), where J P t{v) — S q l' l q| 2 ^( r 7 — ^q) 
is the relevant spectral density for the photon environment. Thus, the spontaneous emission rate we derive within 
the variational theory is dependent upon both the phonon and photon bath correlation functions, C(r) and X(r), 
respectively, and whether this rate varies from that in the absence of the phonon environment depends crucially on 
their respective timescales. 

Consider now the case of no phonon environment, such that we can set C(r) — > 1. Then, we have 

f-t poo 

ri(t) = 2Re / dre iW0T / d ve -" JT J pt { V ), (19) 
Jo Jo 

which, if take the limit of t — > oo, gives the standard Markovian spontaneous emission rate 1^ = 27rJ pt (u;o), dependent 
upon the photon bath spectral density evaluated at the excitonic energy splitting cjo- 

Now, consider the opposite limit in which the exciton-phonon interaction is such that the QD polaron states are 
reached on an extremely short timescale. In this situation, we replace C(r) — > C(oo) = B 2 in Eq. (fT8]) . and we then 
obtain the modified spontaneous emission rate 

r' = B 2 T 1 = 2B 2 7rJ pt (u ), (20) 

again within the Markov approximation. 

We now return to the full expression for T(t) in Eq. (TT8")) to determine the appropriate form for the present situation 
of a coherently-driven driven QD coupled also to an acoustic phonon environment. In this case, we know that the 
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typical timescale for the phonon-bath correlation function to reach the long-time value of B 2 is of the order of a few 



picoseconds [13l. Il4j. For the radiation field correlation function, X(t), we take the standard (3D) spectral density 



J pt (?y) = Arfe 7? / 7? °, where a high-frequency cut-off r) c has been introduced. This gives 

a A 

X{r) = _° , (21) 
(Vc + irp 

which decays to zero on a timescale of roughly l/rj c . Hence, in the limit of large rj c , this correlation function becomes 
extremely short-lived [? ]. We know that in the present case, r\ c must be at least as large as u>o (the ground-state to 
single-exciton state energy splitting), otherwise spontaneous emission would be suppressed. Thus, we can estimate 
i] c > 1.5 x 10 3 ps _1 , for a typical |0) to \X) energy splitting of 1 eV, which leads to a radiation field correlation 
time of the order of femtoseconds at most. On this timescale, the phonon correlation function barely changes, and we 
may replace C(r) by C(0) = 1 in Eq. (|18[) . and we are also now justified in taking the upper limit of integration to 
infinity for timescales of interest. We therefore conclude that, for a typical QD system described by the Hamiltonian 
of Eq. ([5]), the spontaneous emission process is unaltered by the exciton-phonon coupling, and can be described by the 
standard Lindblad form of Eq. (fT5|) with T(t) replaced with Ti = 27rJ pt (u;o). We can thus write the full variational 
frame interaction picture master equation as the sum of two terms 

<I, " V!, ' ) =/C ph (Mi))+ri (a-p v (t)o- + -l{a- + a-,p v (t)}) , (22) 



dt x v 2 

the first describing exciton-phonon interactions within variational theory, and the second giving rise to spontaneous 
emission processes. 



Phonon coupling terms 

Having separated the phonon and photon coupling contributions, we now use the methods described in Ref. [32| to 
find the form of the phonon terms. The variational parameters upon which H\ a , Hit, and H$ all depend are found 
by minimising a free energy bound on the interaction terms. We find /k = g^FlLu^) with 

. (l^tanh(/W2)) 

F K) = t - 2 v-, (23) 

1 - £ tanh(/V2) (l - life coth(^ k /2)) 

and we note that this means the renormalised driving strength, 



Q r = Q exp 



2 



F(w) 2 coth(/3w/2) , (24) 



and R — ^ k w k 1 /k(/k — 2^k) must be simultaneously solved for self-consistently. 

In the Schrodinger picture, pv(t) = exp[—iHst]pv(t) exp[iHst] — i[Hs, pv], the phonon terms can be written 

£ P h (Pv (*))=- 3 S £ ^ ^ ^ ' A = ^ pv W _ pv W A i ^ 

ij u 

- 1 E E S » M > A i (^P v (*) + ^ (t) A ] Ml . (25) 



where {i,j} e {1,2,3} and uj £ {0,±£}. We define A x = a x , A 2 = <r y and A 3 = (1/2) (J + a z ), while A x (0) = 
sin20(|+)<+| - l-X-l), A x (0 =cos2^|-)<+|, A 2 (0) = 0, A a (0=i|-)<+|, 4»(0) = cos 2 |+)(+| + sin 2 9 \-){-\ and 
Az(£) = — sin6>cos6> |— )(+|, defined in terms of the eigenstatcs of H$, satisfying Hs |±) = (l/2)(i? ± £) |±). In all 
cases Ai(uj) — A\{—uS). The angle 6 = (1/2) arctan(il r /e) characterises the tilt of the system eigenstates away from 
the x-axis in the variational frame. Eq. (|25j) contains the quantities = 2Ke[Kij(u>)] and Sij(u>) = Im[l£y(u;)], 

defined in terms of the response functions 



/•OO 

/ AyMe^dr, 
Jo 



Kt^u) = / A ii (r)e^*dr, (26) 
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which themselves depend on the bath correlation functions Ay*(r) = Ti(Bi(T)Bj (0)ps) . Note that in Eq. (|26l) we have 
extended the u ppe r limit of integration to infinity which, for the parameters considered in the main text, is a good 
approximation [28|. We label the bath operators Si = (Sl/2)B X , B 2 = (Sl/2)B y , and B 3 = X)k(#k — /k)(^k + ^k)- The 
bath correlation functions are found to be A u (r) = (^/8)( e ^ T ) + e^M - 2) and A 22 (r) = (Sl 2 r /S)(e^ - e'^), 
with phonon propagator 

' dw ^2 FH 2 G+(T)j (27) 

defined in terms of G±(r) = (n(oj) + l)e~ lWT ± n(w)e JWT , with n(u>) = (e^ w — 1) _1 the occupation number, while 

/•OO 

A 33 (r)=/ dujJ(u)(l-F(uj)) 2 G+(T), (28) 
Jo 

A 32 (r) /°° dw^f («)(1 - F(w))iG_(r), (29) 

2 Jo w 

with A 32 (r) = -A 23 (r), and A i2 (r) = A 2 i(r) = Ai 3 (t) = A 31 (r) = 0. 

Putting everything together, we arrive at the full variational frame Schrodinger picture master equation 

= ~^t ecr * + to r cTx, Pv(t)] +lC ph (pv(t)) +Ti (a-p v (t)<T+ - ^{a + a-,p v (t)} \ , (30) 

as used in the main text. We note that in moving the spontaneous emission terms back into the Schrodinger picture 
the QD operators a± have remained unchanged to be consistent with the approximation leading to Eq. (|11[) . Making 
use of the quantum regression theorem in the standard way, Eq. pop can then be used to calculate the field correlation 
properties of the QD. Slight care must be taken, however, when relating QD expectation values in the variational 
frame to those in the original frame. If \ and xv are the original and variational frame total density operators, 
respectively, we have tr s+B [xcr z ] = tr s+ B[e~ V Xve V &z] = ^s+bIxv^z]- However, for n = x,y, we have tr s+B [xcr„] = 
tr s+s[ e_V Xv® V &n] = Btx s+BlftvPn] where we have made a Born approximation in the variational frame, xv ~ PvPb- 



Pure-dephasing limit 

Though we use Eq. (|2"5|) to numerically calculate the field correlation properties of the QD, the analytical expressions 
resulting from it are somewhat cumbersome. However, in the correct (weak-driving) limit, we find that the phonon 
coupling terms in Eq. (1251) can be well approximated by a simple pure dephasing form. Assuming now that we drive 
the QD on resonance with the polaron shifted transition frequency, v — /Ck^ ^k = 0> anc ^ we drive weakly enough 
such that £1 <C lo Ci then the variational transformation reduces approximately to the full polaron form, and we can 
thus set F(urk) ~ 1 in Eq. (f23| . As such, we find that only the correlation functions An (t) and A 22 (r) survive, and 
the variational master equation [Eq. ([25]) ] reduces to the polaron form given in Ref. [28] . This corresponds to Bloch 
equations of the form a — M ■ a + b where 

f-(T z -T y ) \ 
M = -T y -O r , (31) 

\ o (sir + A) -r, / 

and b = (— k x ,0,0) t , with Bloch vector a = ((a x )t, (& y )t, (c z )t) T . Here, the rates and energy shifts are given by 

T y = 2 7 n(0), (32) 
= (722(n r )+7 2 2(-O r )+2 7u (0)), (33) 
A = 2(S 22 (Sl r )~S 22 (-n r )), (34) 

K = (722 (Sl r ) - 722 (Sir)) ■ (35) 



We are interested in the solutions to these Bloch equations for arbitrary initial conditions as, with the help of the 
regression theorem, this will fully determine the first-order field correlation function and hence the QD emission 
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spectrum. We find that the solutions to the Bloch equations can be written 



(<7z)t 



v — r 



e -(r„+r < )t/2 
e -(r B +r 8 )t/2 



2C 



[K)o(2C cos (Ct) + (r 2 - r„) sin (Ct)) - 2O r (a,) sin (Ct)] , 
[(a 2 )o(2C cos (Ct) - (T, - r„) sin (Ct)) + 2(A + O r )(a J/ ) sin (Ct)] 



(36) 
(37) 
(38) 



where C = ^^(A + ft r ) - (l/4)(r z - T a ) 2 . 

For the appropriate phonon spectral density used in the main text, J p h(w) = au 3 exp[-(u/w c ) 2 ], we find that 
in the regime that ksT < lu c (equivalently (3uj c > 1), as the driving strength becomes small, then 711 (0) becomes 
negligible in comparison to 722(±^r)- Hence, T y — > 0, while T z — > (722(^1-) + 722(— ^r))- Additionally, A <C Q r in 
this regime, and (T z — T y )/( s» Y z j ^J~?tf~- ~ (l/4)rf is very small as well. If we additionally impose Q/3 <C 1, such that 
k/T z — > 0, then we may now approximate the Bloch equation solutions as 



e Tzt (cr x } , 



-r»t/a 



,-r»*/2 



(<T y ) cos (Ct) - -^(cr z ) sin (Ct) 
(ct z )o cos (Ct) + -£-{(Ty)a sin (Ct) 



(39) 
(40) 

(41) 



where C — > \/^r (1/4)T|. Now, if we identify jpd = T z , then these are precisely the solutions we expect from a 
simple pure dephasing master equation of the form 



P = 2~lo*.Pj + —2~( azP<Tz ~ p '- 



(42) 



in the relevant regime of 7po/\/flJ — (l/4)7p D being small. Furthermore, if we wish to ensure that the system tends 
to the correct steady state in the long-time limit, we then need to add a term (iK/A)[a y , {a Zl p}] to the right-hand-side 
of Eq. (j42|) . such that the solution for (a x )t becomes (a x )t = e r - [{cx)o^z + n] — This justifies the forms given 
in the main text for weak driving, and is further confirmed by the agreement we see with the full numerical solution 
of the variational master equation in the appropriate regimes. 



